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I.  INTRODUCTION 


The  concept  of  the  two  stage  free  electron  laser  was  advanced  by 
Elias1  as  well  as  others  about  four  years  ago.  In  the  two  stage  free 
electron  laser,  a  conventional  magnetostatic  wiggler  is  used  to 
generate  radiation  in  the  far  infrared  or  submillimeter  range;  and  this 
radiation  is  used  as  a  pump  to  produce  electromagnetic  radiation  in  the 
visible  or  UV.  The  advantage  of  this  concept  lies  chiefly  in  the  use 
of  low  energy  electrons  (<  10  MeV)  whose  energy  can  be  recovered  very 
efficiently  after  FEL  interaction.  In  our  last  study2,  we  investigated 
several  methods  of  reducing  electromagnetic  pump  optical  cavity  losses. 
The  general  design  considerations  for  a' two-stage  free  electron  laser 
were  given.  A  point  design  for  the  proposed  UCSB  experiment  was  given. 


In  this  study,  we  have  tried  to  validate  the  earlier  point  design 
with  actual  electron  dynamics  calculations.  We  find  the  agreement 
between  the  exact  one-dimensional  code  and  the  point  design  parameters 
to  be  satisfactory.  The  calculations  are  carried  out  for  the  first 
stage  interaction  only.  We  have  also  developed  a  program  to  calculate 
the  transverse  mode  structure  of  the  pump  radiation  in  the  cavity. 
This  will  be  used  in  future  calculations  for  assessing  the  performance 
of  the  second  stage  of  the  FEL.  In  Section  II,  we  have  summarized  the 
results  of  design  and  code  verification.  In  Section  III,  we  discuss 
the  transverse  mode  structure.  Section  IV  contains  a  summary  of  the 
results  so  far  obtained. 


II .  ELEMENTARY  CONSIDERATIONS 


A  conceptual  design  of  the  first  stage  FEL  cavity  was  carried  out 
previously.  This  involved  compromising  on  the  gain,  saturation  flux 
etc.  High  gain  usually  implies  low  saturation  flux  and  vice  versa. 
One  needs  high  gain  to  start  the  first  stage  FEL  quickly,  but  at  the 
same  time  one  also  needs  high  saturation  flux  since  that  will  be  the 
pump  field  for  the  second  stage.  In  this  section  we  shall  first  brief¬ 
ly  go  through  some  of  the  scaling  relations  that  determine  the  wiggler 
parameters.  These  have  been  described  in  detail  before2.  Having 
chosen  a  set  of  design  parameters,  we  shall  describe  the  electro¬ 
dynamics  code  that  actually  solves  for  the  gain  of  the  FEL. 

We  shall  first  assume  that  the  first  stage  wiggler  is  made  out  of 
SmC05  permanent  magnet  as  shown  in  Figure  1.  For  a  magnet  bar  height  g 
*  3  Av/8  and  a  fill  factor  of  unity,  i.e.  no  gaps  between  neighboring 
magnets,  the  normalized  wiggler  vector  potential  on  axis  can  be 
written3 

a  -  1.07  x  lO**4  B  X  exp  [-2*  h/X  ]  (2.1) 

w  r  w  1  wJ 

where  By  is  the  remnant  field  in  gauss  and  other  quantities  are 
defined  in  Figure  1. 

The  photon  beam  width  at  the  entrance  to  the  wiggler  is  denoted  by 
2w  and  for  reasonable  clearance  the  magnet  separation  2h  can  be  taken 
to  be  at  least  a  factor  of  a  larger  than  2w.  The  Rayleigh  range  Zg 
’ is  related  to  the  photon  beam  waist  radius  w0  by 

Z  ■  *  w2/*  (2.2) 

R  O 

The  length  of  the  wiggler  Ly  is  taken  to  be  2  q  times  the  Rayleigh 
range.  Thus  the  beam  radius  at  the  entrance  to  the  wiggler  can  be 
written  as 


w  =*  w  [l+(L  / 2ZD)]1/2 

0  "  R  <2. 3) 

-  [XZr  (Uq2)/*]1'2 

Replacing  h  in  equation  (2.1)  by  aw  and  substituting  (2.3)  we  get, 

1/2 

a  -  1.07  x  10-4  B  X  exp  [-2ira  [X  Zn(l+q2)/ir]  /X  ]  (2.4) 

w  r  w  L  1  R  J  wJ 

The  output  wavelength  X  is  related  to  Xw  and  a,,  through 


X  -  X  (1  +  a  2)/2y2  (2.5) 

w  w 

The  two  convenient  sets  of  independent  parameters  are  Xw  and  Zr  and 
we  shall  plot  the  different  physical  quantities  in  maps  with  Xw  and 
Zr  aa  axes.  Figure  2  shows  a  plot  of  constant  lines.  In 

plotting  the  figure,  we  had  to  assume  a  value  for  y.  The  y  that  can  be 
obtained  for  the  UCSB  experiment  is  12  and  we  have  taken  that  to  be  a 
fixed  value.  Changing  y  would  result  in  similar  plots. 


The  small  signal  gain  g  L 

o  w 


per  pass  is  given  by 


4 


8oLw 


r2  n  B2X2 

o  e  w _ 

me2  (l  +  a2)3/2 
w 


3/2 

X 


X"2 

w 


L 

w 


(Au/u)2 


(2.6) 


where  r  is  the  classical  electron  radius,  n  is  the  electron  density 
o  e 

in  the  beam  and  (Au/o)  is  the  fractional  line  width  and  is  taken  to  be 
~  1/2N  (N  being  the  number  of  wiggler  periods).  In  Ref.  2,  Eq.  (4.6) 
the  factor  (1+a^2)3^2  has  been  mistyped  as  (l+aw2)1/2.  If  we  assume 
that  the  photon  beam  waist  radius  and  electron  beam  radius  are  matched, 
equation  (2.6)  can  be  rewritten  as 


g  L 
°o  w 


-2 

-  1.142  x  10  0  (■ 


aw  ZR 


2 

-)  (' 


udi7) 

v 


(-^)/amp 


(2.7) 
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where  I  is  the  total  electron  beam  current  and  0  3  v/c.  Once  again, 
this  equation  differs  slightly  from  eq.  (4.7)  of  Ref.  2  which  is 
incorrect.  The  saturation  flux  for  a  constant  parameter  wiggler  is 
given  by 


sat 


=  2.67  x  107  ( 


X2 

w 


a2 

w 


ZR 


(2.8) 


For  a  given  set  of  X  and  Z  ,  a  can  be  determined  from  eq.  (2.4). 
Using  eqs.  (2.7)  and  (2.8),  we  can  then  find  the  small  signal  gain  and 
the  saturation  flux.  Figures  3  to  5  show  contours  of  constant  RQ1^> 
$  and  output  wavelength  plots  for  q  *  1  and  q  *  2.  In  all  the  cases 
the  magnet  gap  to  optical  beam  diameter  at  entrance  and  exit  to  the 
wiggler  was  taken  to  be  equal  to  2.  For  the  UCSB  experiment,  the  same 
electron  beam  is  planned  to  be  used  for  the  second  stage  also.  Thus, 
if  we  fix  y  at  12  and  pick  the  first  stage  parameters,  the  second  stage 
output  wavelength  is  predetermined.  These  are  also  indicated  in  Figure 
5.  Now  we  are  in  a  position  to  specify  the  parameters  for  the  proposed 
UCSB  two  stage  experiment.  The  Van  de  Graaff  generator  is  expected  to 
provide  10  maperes  current .  The  bandwidth  requirement  on  the  pump  wave 
to  be  a  coherent  wave  for  the  second  stage  is  that  AX/X  ~  10-i* .  The 
spontaneous  emission  from  the  magnetic  wiggler  is  given  by  5 


dl(a>) 

dfl 


2.94  x  10-8 


hv 


ev 


i  .  N  watts/ev  of  bandwidth 
ma 

/ma/milliradian2/cm 


(2.9) 


For  a  pump  wavelength  of  230u  and  N  ■  5,  the  spontaneous  intensity 
within  the  bandwidth  of  interest  is  found  to  be  ~  6.7  x  10-8  watts/cm2. 
One  would  need  ~  35  exponentiations  to  get  to  saturation  levels  of  108 
to  10  8  w/cm2. 


In  Figure  6  we  have  plotted,  once  again,  lines  of  constant  gain, 
cavity  flux,  output  pump  wavelen^  .h  and  a  in  the  X  -  plane.  For 
Y*12,  the  pump  wavelength  is  limited  to  230  to  460p  if  we  stipulate 
that  the  output  wavelength  is  in  the  visible  (4000  to  8000  A*).  One 


Map  of  Constant  Small  Signal  Gain  Contours 
for  Higgler  Length  Equal  to  Four  Times  the 
Rayleigh  Range 


would  also  need  a  small  signal  gain  of  at  least  5Z  per  pass  (requiring 
~  700  passes  for  saturation).  This  is  important  as  it  directly  impacts 
on  the  pulse  length  necessary  for  a  successful  two  stage  operation. 
From  the  point  of  view  of  wiggler  fabrication,  one  would  want  to  keep 
aw  <  1.  To  achieve  reasonable  second  stage  gain,  we  should  have  the 
cavity  flux  greater  than  ~  108  w/cm2.  The  shaded  area  in  Figure  6 
indicates  the  available  parameter  space  for  operation.  The  actual 
point  that  one  would  pick  in  this  region  depends  on  the  pulse  length  of 
the  electron  beam  available.  If  we  take  the  cavity  length  to  be  ~  3 
meters,  one  round  trip  time  is  20  nsec.  Thus  for  700  passes,  the  pulse 
length  necessary  is  14  nsec.  To  reduce  the  pulse  length,  we  have  to 
pick  an  operating  point  with  higher  small  signal  gain  at  the  expense  of 
decreasing  saturation  flux.  This  would  reduce  the  second  stage  gain 
and  so  we  cannot  arbitrarily  increase  the  small  signal  gain.  We  can, 
however,  locate  the  first  stage  interaction  region  away  from  the  beam 
waist  and  hope  to  operate  the  second  stage  interaction  region  at  the 
beam  waist.  Usually,  the  advantage  gained  by  such  tricks  is  not  over-' 
whelming.  If  we  choose  gQL^  to  be  10Z  per  pass  (for  10  ampere 
current),  then  for  a  3  meter  cavity,  we  would  need  ~  7  usee  long  pulse 
(neglecting  optical  losses).  A  2  outer  cavity  would  reduce  this  to  ~  5 
usee.  Using  Figure  6,  we  can  pick  a  conceptual  design  point  for  the 
operation  of  the  two  stage  FEL.  It  must  be  borne  in  mind  that  the 
boxed  region's  boundary  is  somewhat  fussy  and  arbitrary.  The  require¬ 
ment  on  the  output  wavelength  can  be  considerably  relaxed  if  one  uses 
two  separate  electron  beam  sources  for  the  first  and  second  stage.  The 
limit  on  the  minimum  gain/amp  can  be  relaxed  if  the  current  can  be 
increased  to  say  100  amperes.  For  the  conceptual  design  of  the  two 
stage  FEL,  we  have  chosen  the  nominal  parameter  values  given  in  Table  I 
for  the  UCSB  experiment. 

For  a  10  usee  long  pulse,  the  fluence  at  the  beam  waist  is 
~6  kJ/cm2 .  An  area  expansion  of  ~  100  would  bring  the  mirror  fluence 
to  60  J/cm2  and  energy  absorbed  with  0.2Z  absorption  would  be  ~  120 
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TABLE  I 

WIGGLER  PARAMETERS  FOR  THE  TWO-STAGE  FEL  EXPERIMENT 


1 

■ 

12 

I 

- 

10  amp 

A 

. 

4  cm 

w 

ZR 

- 

6  cm 

Si 

m 

24  cm 

a 

m 

0.855 

w 

*0^ 

m 

7Z  per  pas  8 

- 

5.8  x  10®  w/cm2 

at 

■ 

2.14  mm 

0 

*P 

- 

240  u 

A 

. 

419  nm 

Number  of  periods  *  6 


-12 


mj/cm2,  which  should  be  quite  safe  in  terms  of  mirror  damage.  If  Lc 
is  the  length  of  the  optical  cavity  and  Rj  and  R2  are  the  mirror  radii, 
we  can  define  quantities  g^  and  g2  by 

gl  2  *  1— L/R^  2  (2.10)) 

»  » 

For  stability6,  we  require  that  0  <  gjg2  <1.  We  shall  consider  the 
symmetric  case  with  gj  *  g2  ■  g.  Since  we  want  the  beam  to  expand 
rapidly,  we  shall  choose  a  design  where  the  cavity  is  nearly  concentric 
(L  *  2R) . 

Let  L/R  ■  2  -  e.  Then 


g  *  e  -  1. 


The  radius  of  the  beam  at  the  waist  is  given  by6 


(  XLc  }1/2  ^8l82  (l  "  *1*2)]1/4 

*  [Si  +  82  "28l  82 l1 /2 


XL  1/2  l'*» 

.( - i)  ( — £_) 

v  2*  ’  v2-e  ' 


Equation  (2.11)  can  be  rewritten  as 


2  L 

*■  R  c  ,  e  x 

*  wo  "  ~1 - IT  (T^) 

R 


1/2 


(2.11) 


(2.12) 


2 

Since  Zg  is  defined  to  be  ir  wq/X,  we  have, 

r  ,1/2 

qc  [c/(2-e)j  -  2 


(2.13) 


-13- 


where  q  *  L  /Z„. 
c  c  R 

The  spot  size  on  the  mirrors  is  given  by 

XL  1/2 
W1  -  (-^) 


«2 

r  _  i 

L  gi(i-gig2)  J 


(2.14) 


(2.15) 

(2.16) 

(2.17) 


Knowing  q,,  and  e,  we  can  obtain  the  cavity  length  and  mirror  radius 
of  curvature.  The  diameter  of  the  mirror  can  be  calculated  by  fixing 
the  diffraction  loss.  Li7  has  numerically  calculated  the  diffraction 
losses  for  different  |g!  values,  and  Fresnel  numbers.  In  general,  if 
Igl  <  0.99,  one  would  need  a  Fresnel  number  of  15  or  larger  to  keep  the 
diffraction  losses  per  pass  to  <  10"4 .  In  Table  II  two  different 
cavity  designs  are  given  for  an  area  expansion  factor  of  100  and  200 
for  a  Fresnel  number  of  15  (■  a2/\Lc,  where  a  is  the  radius  of  the 
mirror). 

The  choice  between  the  first  and  second  design  would  depend  on 
other  factors  such  as  space  to  bend  the  electron  beam  into  and  out  of 
the  wiggler,  electron  beam  pulse  length,  etc. 
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TABLE  II 

OPTICAL  CAVITY  PARAMETERS 


Beam  Area  Expansion 

100 

200 

e 

0.02 

0.01 

qc 

19.90 

28.21 

Cavity  Length  L 

c 

119.4  cm 

169.3  cm 

Radius  of  Curvature 

60.3  cm 

85.1  cm 

Mirror  radius  a 

6.56  cm 

7.81  cm 

Pump  wavelength 

240  pm 

240  pm 

Cavity  Round  Trip 

~  8  nsec 

~  11.3  nsec 

Having  chosen  the  design  values  for  the  wiggler,  the  single 
particle  electron  dynamics  code  was  run  to  verify  the  gain  and 
saturation  characteristics.  The  electron  dynamics  code  solves  the 
following  two  equations: 

1_  (Y|)  .  [g  +  |  x  (?  +  $  )]  (2.18) 

dt  me  L  r  r  w  J 


(2.19) 


where  E  and  B  are  the  electric  and  magnetic  fields  of  the  radiation 
r  r 

and  is  the  wiggler  field.  The  equations  are  solved  for  different 
initial  phases  of  the  electron  with  respect  to  the  radiation  field  and 
averaged  over  the  initial  phase  angle. 


The  results  are  shown  in  Figure  7,  in  which  the  gain  is  plotted  as 
a  function  of  the  cavity  flux  for  different  energy  off-sets  from  the 
resonant  value.  Since  the  energy  spread  in  the  incoming  electron  beam 
is  expected  to  be  less  than  0.  IX ,  it  is  a  good  approximation  to  con¬ 
sider  the  electrons  to  be  monoenergetic .  He  find  the  agreement  between 
the  exact  calculation  and  the  chosen  design  values  to  be  reasonable. 

For  a  fixed  incoming  beam  energy  and  wiggler  configuration,  the 
wavelength  at  which  maximum  gain  occurs  changes  with  increasing  flux. 
If  there  is  no  wavelength  discrimination  done  by  the  optical  cavity, 
the  output  wavelength  will  progressively  shift  towards  longer  wave¬ 
length. 

The  one  dimensional  calculation  described  above  is  correct  for 
only  incoming  plane  waves.  It  is  of  interest  to  study  the  transverse 
mode  structure  under  steady  state  conditions  for  a  given  cavity  con¬ 
figuration  and  this  is  treated  in  the  next  chapter. 
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III.  OPTICAL  MODE  STRUCTURE 

The  equations  that  describe  the  dynamics  of  the  free  electron 
laser  have  been  given  by  many  authors  before8 .  For  the  sake  of 
consistency  in  the  notation,  we  follow  the  work  of  Kroll,  Morton  and 
Rosenbluth9  and  give  their  derivation  below.  We  assume  the  wiggler 
magnetic  field  and  optical  field  to  be  derived  from  the  vector 
potentials 


A  *  /2~(mc2/e)a  (z)  $  sin  JZ  k  (z,)dz,  (3.1) 

w  w  1  w  i  i 

o 

A  ■  /J(mc2/e)a  (z)  Q  sin  /*  k  (z.  )dz,  -  u»  t  (3.2) 

s  s  1  s  1  1  s 

o 

where  the  dimensionless  vector  potentials  a  *  eA/mc2  have  been 

introduced.  We  assume  that  the  fractional  change  in  a  ,  a  ,  k  and  k 

w  8  w  s 

is  small  over  a  distance  of  one  period  of  the  wiggler.  Maxwell's  wave 
equation  for  the  evolution  of  the  electromagnetic  vector  potential  is 
given  by 


Writing  kfl(z) 
obtain 


k  +  6k  and 
so  s 


3t2 


substituting  Eq. 


1  ’ 

(3.2) 


in  Eq. 


(3.3) 

(3.3)  we 


2 

/T  (3£— ) 

6 


2m  2^ 

8  A  _ s 

-  ok  a  e,  - 

c  s  s  *  c 


(3.4) 


where  vectors  ^  and  are  defined  by 


( 
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6k  A 


s  a 


where  <  >  indicates  time  average. 
For  a  single  electron  we  have 


4n  * 

c  <Ji-  el‘> 


J1  .  -  ev^  ^  6[x-XQ(t)  ]6[y-yQ(t)  ]6[z-*o(t)  ] 

where 


♦ 


-► 

-eA 

Ymc 


+ 

-eA 
_ _ w 

ymc 


for  A  «  A  .  Hence 
s  w 

*►  a  e^A^tz) 

(Ji‘  *1  )single  particle  “  “y5f -  tc°8  ,('“C08  (*"2  /\<zi>dzi>] 

o 

x  6[x-Xo(t)]6[y-yo(t)]6[z-*o(t)] 

where 

♦  "  /*[*<*'>  ♦  k „(z')]  dz'  -  to  t 

0  W  8  8 


(3.5) 

(3.6) 

(3.7) 

(3.8) 

(3.9) 

(3.10) 

(3.11) 
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and  x  ,y  and  z  are  Che  coordinates  of  Che  electron  at  time  C.  Time 
o  o  o 

averaging  Co  remove  high  frequency  variations,  summing  over  all  the 
electrons  and  averaging  over  the  beam  cross-sections,  we  obtain 


’-4 ^-(527L).  (3.12) 


C  j  e  ) 

'  1*  1  time  average  me  V  Y 

where  n@  is  the  electron  density  and  ( - )  means  average  over  initial  i|/ 

and  energy.  Substitution  of  Eq.  (3.12)  with  Eq.  (3.7)  yields 


s 


(C08l)< 

Y 


)• 


(3.13) 


Again  from  Eq.  (3.4),  we  find 

2w 

- 2-  A'  -  —  <  J..  e->  ,  (3.14) 

C  8  C  J.  Z 


and  proceeding  in  a  similar  manner,  we  find 


u>2 

a- _ 2 — 

a  2cu 

8 

Prom  Eq.  (3.15),  we  easily  obtain, 


dE 
_ 8 

dz 


2ir  n  e  a 
e  w 


(3.15) 


(3.16) 


where  Ea  is  the  electric  field  strength  of  the  electromagnetic  wave 
(in.e.a.u.)  and  g  is  the  gain  per  cm.  We  therefore  obtain, 


g  * 


2ir  V  /  sinji  \ 

10  Ea  {  Y  / 


(3.17) 
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where  J  is  the  electron  beam  current  density  in  amperes/cm2, 
also  easy  to  see  that 


It  is 


6k  ■  g  (cos  iji/y)  /  (sin  (3.18) 

s 

We  write  E  *  A  u(x,y)  where  u(x,y)  is  normalized  such  that 
so  J 

Jj  u(x,y)  u*(x,y)  dx  dy  *  1.  From  the  relation 

//  EgE£  dx  dy  »  Pg  *  Cavity  Power  (3.19) 


we  find 


4ir  P  x  107  \ 

- 1 -  )  (3.20) 

where  Ps  is  in  watts.  Equations  (3.17)  and  (3.18)  with  the  help  of 
the  above  can  be  combined  to  give 

•  ai  _  2x  /  c 
8  ~  1  s  10  (  4irP  x  10; 

x  |^(sin  i j>/y)  -  i  (cos  i>/y) 

Equation  (3.21)  will  be  the  starting 
resonator  mode  structure. 

Various  techniques  exist  for  finding  the  modes  of  a  bare 
resonator.  One  of  the  first,  and  most  successful  of  these  is  the 
Fox-Li10  technique,  where  an  arbitrary  amplitude  fluctuation  is 
propagated  numerically  between  resonator  mirrors  until  an  eigenmode  is 
obtained.  The  usefulness  of  this  technique  depends  strongly  on 
obtaining  a  good  initial  guess  for  the  mode. 


1/2  a  (z) J(x,y) 
w _ 


] 


(3.21) 


point  for  the  analysis  on  the 


Analytic  techniques  also  exist  for  obtaining  the  modes  of  finite- 
mirror  stable  resonators.  For  the  infinite  symmetric  mirror  case, 
the  solutions  are  those  of  the  integral  equation 

u  (x,y)  ■  //  dx'  dy'  u(x',y')  exp  ^ik  [(x-x')2 

+(y-y’)2]/2L  -  ib  (x2  +  y2)  -  ib  (x'2  +  y'2)^  (3.22) 


where 


b  ■  k/2R,  R  is  the  radius  of  curvature  of  the  mirror  and  L  is  the 
separation  between  the  mirrors.  If  R  *  L  (confocal  resonator)  then  the 
equation  becomes 


u  (r)  »  /  dr'  u(r')  exp  [-ik(r .r ' )/l] 


(3.23) 


where  the  integral  extends  over  the  aperture  of  the  cavity.  The 
solutions  of  Eq.  (3.23)  are  the  functions  which  are  their  own  fourier 
transform.  The  solutions  are  the  well-known  Hermite-Gaussian  functions 
for  the  infinite  aperture.  For  the  case  of  finite  mirrors  the 
solutions  are  the  so  called  prolate  spheroidal  wave  functions.  Away 
from  the  mirror  edges  (>  /XL) ,  they  reduce  to  the  Hermite-Gauss 
functions.  These  solutions  were  first  obtained  by  Boyd  and  Gordon11. 

The  problem  at  hand  is  to  find  the  modes  of  a  loaded  resonator 
with  a  gain  profile  of  the  form  given  by  Eq.  (3.21).  The  technique 
that  we  use  is  that  of  the  degenerate  perturbation  theory.  The 
degenerate  perturbation  theory  treatment  is  necessary  since  the  eigen¬ 
values  of  the  lowest  modes  can  differ  from  each  other  and  the  fundamen¬ 
tal  mode  by  less  than  that  of  the  matrix  element  of  the  perturbation12. 
Basically,  we  want  to  take  linear  combinations  of  the  lowest  order 
modes  via  matrix  methods  to  find  a  prototype  "unperturbed  state.”  Once 
this  is  done,  the  ordinary  perturbation  theory  can  be  implemented.  The 


( 
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reason  for  following  Che  perCurbation  theoretical  approach  is  that  the 
integrated  single-pass  gain  of  the  free  electron  laser  is  expected  to  be 
less  than  1,  typically  0.2  to  0.5.  In  cases  when  the  validity  of  the 
perturbation  theory  is  questionable,  the  results  of  this  analysis  can  be 
used  as  the  initial  guess  for  Fox-Li  analysis. 

The  geometry  of  the  problem  is  sketched  in  Figure  8.  We  propagate  a 
ray  of  light  from  a  point  r  in  the  plane  z  ■  0  through  the  gain  medium  to 
the  right  mirror,  then  back  to  the  left  mirror  and  finally  back  to  the 
point  r*  in  the  mid-plane  z  ■  0.  It  is  assumed  that  the  electrons  move 
from  left  to  right .  It  must  be  borne  in  mind  that  the  FEL  gain  is 
unidirectional.  The  integral  equation  for  the  loaded  system  involves  the 
propagator  between  the  points  r  and  r'  in  the  figure.  The  equation  is  of 
the  form 


y  u  (x,y)  -  /  R  (r  ,r 1 )  u  (x',y')  dx'  dy' 
mm  m 


(3.24) 


where 


K  -  K(0)  +  -  K<0)  exp  {  J*  ( g  -  i  6ks>  ds} 


(3.25) 


(°)  . 

over  the  path  shown,  y  is  the  eigenvalue  and  K  is  the  kernel  of  the 

m 

propagator  for  the  bare  resonator  given  in  Eq.  (3.23).  If  /  (g-i6k) 
ds  «  1, 

(1)  (0) 

then  K  *  K  J+  (g-i  Sk^)  ds .  Actually,  the  approximation  is 


better  than  this,  since  the  gain  g  for  the  variable  parameter  FEL 

(°) 

is  linear  rather  than  exponential.  K  represents  the  effects  of  the 

ordinary  Fresnel  propagator  and  mirror  curvature.  To  simplify  the 

algebra  and  the  notation  in  what  follows,  we  shall  use  matrix  form  and 

use  K  interchangeably  with  the  integral  operator  and  the  kernel  R(r,r'). 

(0) 

The  eigenvalues  and  the  eigenfunctions  of  the  operator  K  are  known  (as 
given  by  Boyd  and  Gordon)  for  a  given  set  of  cavity  parameters  and  we 
take  these  to  be 
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(3.26) 


0) 
u 

n 


where  u  are  the  prolate-spheroidal  wave  functions  as  mentioned  before 

and  Yn  are  the  corresponding  eigenvalues.  The  superscript  (°)  refers 

to  the  order  of  the  solution  and  the  subscript  'n'  refers  to  the  n^ 
eigenvalue  or  eigenfunction.  To  apply  the  degenerate  perturbation 

theory,  we  shall  assume  s  of  these  eigenvalues  to  be  degenerate  or  near 

(°) 

degenerate.  The  choice  of  the  eigenfunctions  ur  for  the  degenerate 

level  is  usually  arbitrary,  but  ceases  to  be  so  if  we  subject  the 
eigenfunctions  to  the  requirement  that  the  change  in  them  under  the 
action  of  small  applied  perturbation  be  small.  Let  the  correct  choice  of 

s  (°)  (°) 

the  new  basis  of  eigenfunctions  be  of  the  form  E  C  u 

m»l  m  m 

We  write 


(°)  (l) 

K  -  K  ’  +  K  ♦  ., 


_  (°)  (!)  ^ 
u  ■  u  +  u  +  •  .  . 
Hi  n  n 


(°)  O) 

v  ■  v  +  y  ▼  • . . 

m  m  m 


(3.27) 


Straightforward  application  of  the  near  degenerate  perturbation  theory 
yields 


l 

n 


It  „(1)  -CT-r  )(0) 


6 

mn 


0 


(3.28) 
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where 


(!) 


am 


-  //  un(0)  (x,y)  K(1)  (x,y ;  x\y')  um(0)(x\y') 
dx  dy  dx'  dy' 


(3.29) 


For  this  system  of  homogeneous  linear  equations  for  the  quantities 

(°) 

C  to  have  nontrivial  solutions,  we  require  that  the  secular  equation 


-  (r-Y  (0))  I 


(3.30) 


be  satisfied. 


Substituting  in  turn  the  roots  of  Eq. 

and  solving,  we  find  the  coefficients 
eigne function  in  the  zeroth  approximation. 


(3.30) 


in  the  system  (3.28) 
and  so  determine  the 


The  first-order  correction  to  this  eigenfunction  in  the  orthogonal 
space  is  obtained  just  as  in  regular  perturbation  theory.  The 
.  .  (!) 

coefficients  ,  where  k  *  m,n,...  are  obtained  as 


n  n 
Y  -  Y 


(3.31) 


The  complete  solution  correct  to  first  order  is  given  as 


un(x’y)  "nil  Cn(0)  un<0> 


(°)  O) 

•  a  C  K. 

*  k-5+i  A  -777^-  \ 


(3.32) 


For  numerical  calculations,  the  series  in  the  second-term  is  truncated  at 
a  convenient  point  determined  by  the  accuracy  of  the  procedure.  This 
then  completes  formal  solution  to  the  problem.  For  the  FEL,  the 
matrix  element  kJ  >  can  be  written  as 


Kmn  “  AL  ^  dx  dy  dx'dy'  ^(x,y).  u^^x’.y’) 

*  exp  (i2b  [<x'-x3)2  ♦  <y*-y3>2]  -  ia  [^2  +  ^2] 

+  ib  [(*3-x2)2  ♦  (y3-y2)2]  -  ia  [^2  *  ^2] 

*  i2b  [<*2-x)2  +  <y2-y)2J}  s  (x,y;x',y-) 


(3.33) 


where 


*2  "  (3x+x')/4 
x3  *  (*+3x')/4 
«  ■  */XR 

9 

the  function  S  (x,y;  x’y‘>  is  given  by 


y2  -  (3y+y 1 )/4 
*  (y*3y')/4 
b  -  ir/XL. 


•  r.x'.y')  -  /°  dz  [u  jx-x'^-Kw'^,  1/2  I 

-in  41,2  1 

*  X 

*  {  K'2“  )»«  (,.iSk>)  j 


(3.34) 
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where 


*i(z)  “ x’ +  Hhr^  z ;  yi(z)  3  yl  +  (_z5?)  z 

x2(z)  3  X  ♦  (  ~2^~)  Z  ;  y2(z)  3  y  ♦  *•  (3.35) 

The  electron  and  photon  beam  radii  at  the  waist  are  typically  of  the 

order  of  1  ■  while  the  mirror-separation  may  be  ~  1  meter  and  the 

wiggler  length  1^  may  be  ~  0.2  meters.  We  therefore  expect  z/2L  to 

have  a  maximum  value  <  0.25  and  one  can  therefore  expand  g(x1,y1)  around 

g(x',y')  in  a  Taylor  series.  Similar  expansion  holds  good  for  g(x2,y2). 

The  z-integration  can  be  trivially  carried  out  to  the  desired  order  in 

(1^,/L).  The  integral  in  Eq.  (3.33)  is  a  four-fold  integral  which,  in 

general,  cannot  be  reduced  to  simpler  integrals.  The  numerical 

(1 ) 

evaluation  of  K  is  probably  most  conveniently  carried  out  by  Monte- 

mn 

Carlo  methods.  The  algebra  can  be  considerably  simplified  if  the  gain 
can  be  written  as  a  product  of  separable  functions  in  x  and  y.  The 
current  density  J(x,y)  in  Eq.  (3.21)  will  be  taken  as  Gaussian  while  for 
u(x,y)  in  the  same  equation,  we  can  take  the  Gaussian  form  as  the  first 
guess  for  the  solution.  This  will  enable  us  to  break  the  four-fold 
integral  in  Eq.  (3.33)  to  the  product  of  two  double  integrals  which  can 
be  conveniently  carried  out  in  a  computer. 

The  numerical  calculations  that  have  been  carried  out  incorporate 
further  simplifications.  The  zeroth  order  eigenfunctions  were  taken  to 
be  Hermite-Gaussian  functions  rather  than  the  prolate  spheroidal 
functions  since  the  two  differ  only  at  the  edges  of  the  mirror  and  the 
former  is  easily  calculated  numerically.  It  must  be  emphasized  that  this 
simplification  is  not  essential  to  the  method  outlined  here.  Secondly, 
the  zeroth  order  eigenvalues  were  calculated  approximately  by  the 
following  means.  We  first  calculate  quantities  B2  defined  by 
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(3.36) 


(J2 

m 


/  in  IH  (JS.)  1 2 exp  (-2xV»,2]  dx 
-d/2  ^  Wl  _ 

/"  lHn  (“^p}  !Z  «p  [-2x2/Wl2]  dx 


where  w,  is  Che  beam  radius  at  Che  mirror.  The  quancicy  02  62  gives  Che 
1  in  n 

f  race  ion  of  power  in  Che  TEM  modes  Chac  is  inCercepCed  by  Che  mirror. 

mn 

The  e  ige  n  f  unc  C  i  on  H  ( — — *  1  H  )  exp  [-  (x2  +  y2)/w  21 

*  m  w  J  nen1-  w  J  1  o  J 

Ch  0  0  (°)  . 

represenCs  some  k  eigenfunction  of  Che  operacor  K  approximately. 

Hence,  we  expect  y,  3  6  6  .  For  Che  numerical  example,  we  cake  L  *  1.2 

lc  o  n 

meters,  L*  *24  cm,  and  R  *  60.3  cm.  In  Figure  9,  we  have  plotted  Che 
intensity  profiles  on  the  mirror  for  different  cavicy  power  levels.  It 
is  interesting  co  see  Chat  Che  four  lobed  distribution,  indicative  of 
TEMqi  and  TEMjq  modes,  present  at  small  power  levels  gives  way  to  a 
central  peaked  intensity  distribution  at  high  power  levels.  If  we  had 
used  circular  mirrors  and  cylindrical  coordinates  instead  of  square 
mirrors  and  cartesian  coordinates,  the  four-lobed  distribution  would  be 
replaced  by  an  annular  (doughnut  shaped)  intensity  distribution. 

For  the  calculation  of  the  performance  of  the  second  stage  of  the 
two-stage  FEL,  it  is  also  of  interest  to  know  the  intensity  distribution 
at  different  zr locations  in  the  wiggler.  These  are  shown  in  Figures  10 
to  16  at  a  cavity  power  level  of  1.44  x  107  watts. 
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Figure  10.  Intensity  Distribution  of  z  ■  -12  cm  for  a  Cavity  Power 
Level  of  1.44  x  10?  watts.  Other  parameters  are  the 
same  as  in  Figure  10.  z  ■  0  is  the  center  of  the  wlggler. 


IV.  SUMMARY 


We-4iau«  verified  Che  simple  conceptual  design  parameters  of  the 
two  stage  FEL  with  detailed  calculations  of  electron  dynamics.  The 
agreement  between  the  two  is  quite  reasonable.  These  calculations 
assume  a  plane  wave  propagating  through  the  wiggler.  For  a  given 
stable  cavity  configuration,  the  actual  transverse  intensity  profiles 
were  then  calculated  using  perturbation  techniques.  The  theory  used 
should  be  accurate  for  single  pass  gains  less  than  50Z  and  for  cavity 
lengths  long  compared  to  the  wiggler  length. 

In  the  continuation  of  this  program,  the  results  of  the  present 
calculations  will  be  used  to  determine  the  performance  of  the  second 
stage  of  the  free  electron-  laser  that  will  produce  output  in  the 
visible  portion  of  the  electromagnetic  spectrum. 
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